Influence of geography on language competition 
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Competition between languages or cultural traits diffusing in the same geographical area is studied 
combining the language competition model of Abrams and Strogatz and a human dispersal model on 
an inhomogeneous substrate. Also, the effect of population growth is discussed, ft is shown through 
numerical experiments that the final configuration of the surviving language can be strongly affected 
by geographical and historical factors. These factors are not related to the dynamics of culture 
transmission, but rather to initial population distributions as well as geographical boundaries and 
inhomogeneities, which modulate the diffusion process. 
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I. INTRODUCTION 



pected effects. Results are summarized in Sec. IVIII 
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Currently, statistical mechanics and stochastic models 
are employed to study a wide range of topics, not only 
in physical sciences, but also in interdisciplinary applica- 
tions in biology as well as in social and historical sciences. 
Examples include applications to financial time series [l| , 
population dynamics [21, and archeology Q- Recently, 
also various problems in linguistics have been approached 
using methods imported from the theory of complex sys- 
tems and statistical mechanics; see Refs. [J, a B fo'" ^^ 
overview. 

In the evolution and dispersal of biological species the 
importance of geography is well known JJ . The goal of 
the present paper is to study how purely geographical and 
historical constraints can affect the dynamics of different 
languages or cultural traits competing in a region. Pre- 
viously the influence of geography on language dynamics 
has been investigated e.g. in Refs. fa. M. [lOl. \V\\ . 

We start from the Abrams-Strogatz (AS) model [IJI of 
two fixed competing languages. The term "fixed" refers 
to the fact that the evolution of language is neglected 
on the time scale considered, so that the model is for- 
mally similar to a model of population dynamics of two 
biological species. By "competing" one means that at 
any time speakers can switch to the other language, as a 
consequence of the interaction between speakers of lan- 
guage 1 and 2. In order to take into account population 
growth, dispersal, and the effect of geographical inhomo- 
geneities, we introduce in Sec. [ll] a more general model. 
This model is then applied to some (idealized) examples 
concerning the influence of initial conditions (Sec. IIIip . 
boundary conditions (Sec. lIV[) . and geographical barriers 
(Sees. IVl and rvT]) . These examples show how extending 
a 0-dimensional (i.e. homogeneous) model of a culture 
transmission to physical space gives rise to new, unex- 



II. THE MODEL 

In this section we formulate a model of culture disper- 
sal and competition. We are interested in a description 
regarding a historical context [20|, in which two popula- 
tions using different languages (or with different cultural 
traits) 1 and 2 disperse across a region. Simultaneously, 
the competition between the languages takes place and 
possibly the populations also grow. As a starting model 
we assume the AS language competition model. A dif- 
fusion term, describing population dispersal, and an ad- 
vection term, taking into account geographical inhomo- 
geneities, are then added; population growth is taken into 
account through a logistic term. 



A. The Abrams-Strogatz model 

The AS model was first introduced in Ref. [I^l in order 
to describe the time evolution of the population size of 
various endangered languages. The model can be formu- 
lated through the following equations, 
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Here Ni{t) {i — 1,2) represents the fraction of speakers 
of language i. The quantity fci = Si/t is the rate con- 
stant for the switch of a speaker of language 2 to language 
1, and vice versa for ^2 = S2/T. The dimensionless pa- 
rameter Si, referred to as language status, represents the 
attractiveness of language i; it can result from a combi- 
nation of factors, such as, e.g., the language prestige and 
usefulness. Here we follow the normalization convention 
S1 + S2 = 1. As a consequence si/t + S2/t= 1/t, i.e., r 
represents the time scale. The coefficient a varies slightly 



around the value a — 1.3 for different languages, as found 
in Ref. [l3|; in this article the value a — 1.3 is assumed. 

In the analogy with a population dynamics model it 
should be noticed that the reaction term R{Ni,N2) in 
Eqs. H]) contains a positive and a negative contribution, 
depending on both populations Ni and N2 , which repre- 
sent an advantage and disadvantage due to the encounter 
with an individual of the other "species" . In other words, 
speakers of population 1 and 2 behave symmetrically to 
each other as prey and predator at the same time [2]. 
While this situation is not usual in biology, it can be jus- 
tified for the interaction between two cultural traits [ij] . 

The analysis shows that the AS model has one unstable 
and two stable equilibrium points. The latter ones corre- 
spond to one of the languages surviving and the other one 
disappearing. Which final state will be reached depends 
on the initial populations Ni{tQ), the status parameters 
Si, as well as on the value of a. The critical values of 
parameters defining the unstable equilibrium point can 
be obtained from Eqs. H]) setting the rate term R equal 
to zero. 
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When the ratio NI/N2 is larger than the right hand 
side of condition ((21) at some time t ^ t' , then 
R{Ni{t), N2{t)) > at any later time t > t' and iV2 ^ 
for t^oo, while Ni->N' = Ni{t')+N2it'). The opposite 
takes place if N^ /N2 is smaller than the right hand side. 



B. Generalized model 

Population and culture spreading may be affected by a 
wide range of geographical factors, due to physical barri- 
ers such as water boundaries and mountains or e.g. geo- 
physical features such as type of ground and spatial dis- 
tribution of resources [7|. While the underlying mech- 
anisms determining the influence of such geographical 
factors are in general complex, in a first approximation 
their overall effect can be described statistically. In fact, 
dispersal of human populations in an geographical envi- 
ronment recalls the diffusion of Brownian particles modu- 
lated by an external field or an inhomogeneous substrate. 
For example, human dispersal in neolithic Europe 13 1 
and during the early colonization of South- America 14 1 
has been studied using advection-diffusion equations with 
a logistic term taking into account population growth 
(i.e., employing the so-called Fisher equation). 

In order to extend the AS model to take into account 
the geographical inhomogeneities, we merge it with the 
two-dimensional geographical model of human dispersal 
and growth proposed in Refs. [IJ, ^^ . The corresponding 



evolution equations read 
dfi 



dt 



dh 
dt 
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with the reaction term given by 

R{h,h)=k{s,f,-f2-S2f2^h) 



(3) 



(4) 



The quantity fi — fi{x,y,t) represents the population 
density of speakers of language i, while the constant k 
in Eq. ([4]) is an effective rate constant and Si still rep- 
resents the status of language i. In Eqs. ^ population 
movement is described by the advection term containing 
the external force field r(a;, y) — {Fx{x, y),Fy{x, y)) and 
by the diffusion term with the diffusion coefficients D = 
D{x, y, ). In general F and D are related; their explicit 
form depends on the problem considered. One possibil- 
ity is to set F{x,y) = and describe the inhomogeneous 
character of the substrate through a space-dependent dif- 
fusion coefficient D{x,y) [iJillJI- ^^ the model systems 
studied below we ascribe the inhomogeneous character 
of dispersal to the external force, F = F(x,y), while D is 
kept constant. For illustrative purposes and in analogy 
with Brownian motion, the force field is expressed as the 
gradient of a potential, F(a;, y) = — VC/(x, y). The logistic 
terms with Malthus rate a and carrying capacity K in 
Eqs. ([3]) take into account the population growth. 

According to Eqs. ([3]), populations 1 and 2 disperse in- 
dependently, whereas the densities /i and /2 are coupled 
only through the reaction term R and through the logistic 
terms, which introduce a negative competitive coupling 
proportional to — /1/2. Also, it should be noticed that 
Eqs. ([3]) describe two populations with identical disper- 
sal and growth properties, differing only in the way how 
the respective language interact, according to the term R 
given by Eq. Q. As a consequence, the total population 
density f = fi + f2 follows a diffusion-advection-growth 
process without culture transmission, obtained by sum- 
ming Eqs. ([2]), 



1^ = -V • [F/] + V(I?V/) +af(l-^ 



(5) 



In order to solve Eqs. ([3]) numerically, one can ap- 
proximate the derivatives through the corresponding fi- 
nite differences, replacing the problem in the continuous 
time and space variables (i, x, y) with the one on a lat- 
tice {k,m,n) defined by the discrete variables tk = k6t, 
Xjn = mSx, and j/„ = n6y, respectively, where fc, to, 
n are integers, while St, Sx, Sy are the lattice steps. By 
using the explicit Euler integration scheme |15j . for con- 
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FIG. 1: Comparison of the evolution of population densities fi{x,y,t) and f2{x,y,t) (columns 1 and 2) for two languages with 
status si = 1— S2 = 0.55 for different widths of the initial distribution f2{x, y, 0). The initial distributions fi{x, y, 0) are given by 
Eq. ^. Example A: a more localized initial distribution f2{x,y,0) with ai = 1.75. Example B: a more spread f2{x,y,0) with 
(7i = 3. All other parameters are the same, see text for details. 



stant D one obtains from ^ the finite-difference equa- 
tions 



St 



= ±R{fuh)- 



6x Sy 



-D 



S^ f S'^f 
5x^ Sy^ 



+ a/. (1-;^) , * = 1,2. (6) 



Here the term -\-R corresponds to i = 1 and —R to i = 
2. The time difference operator Sk, when apphed to a 
generic function 4>{k, m, n), gives 

5k4>(k, m, n) = (j){k + 1, ra, n) — (t){k, m, n) . 

As for the x-difference operators, 5m is defined as 

Sm (t>ik, m, n) — [(t>{k, m+l,n) — (t>{k, ?7i— 1, n)]/2 , 

while the second difference operator S^ as 

S^(f>{k,m, n) — (f){k, ni+l, n) + (j){k, ni— 1, n) — 2(j){k, m, n) ; 



analogous definitions apply to the y-difference operators 
5n and 5^ . If the finite-difference equations (O are used 
to solve numerically the continuous problem, suitable 
constraints have to be imposed on the values of the steps 
5t, 6x, Sy in order to limit the numerical integration er- 
ror [l5|. 

One can also consider Eqs. (O as a general lattice 
model of dispersal, growth, and cultural interaction of 
two populations, with no reference to the continuous 
equations ([3]). Such a discrete model is employed in the 
examples presented in the forthcoming sections, even if 
for convenience we will refer also to the corresponding 
continuous limit represented by Eqs. ([3]). 



III. INFLUENCE OF INITIAL CONDITIONS 

The minimal spatial version of the AS model is ob- 
tained when taking in Eqs. ^ the diffusion coefficient 
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FIG. 2: Time evolution of the total population N\{t) = 1 — 
N2{t) with a status si = 0.55 = 1 — S2 for the examples A 
(continuous line) and B (dotted line) of Fig. [l] The critical 
fraction A''i — 0.3387 given by the AS model for the survival 
of language 1 and the line corresponding to A^i — N2 are also 
drawn. 



D equal to a constant, and neglecting the growth and 
advection terms. Shxiilarly to the original AS model, it 
admits two stable equilibrium solutions, corresponding 
to one of the two languages surviving and the other one 
disappearing, as can be shown by linear stability analy- 
sis. Depending on the initial conditions, there exist also 
unstable equilibrium solutions. In the presence of a local 
spatial noise the unstable solutions are washed out by the 
random fluctuations [10]. However, the minimal spatial 
model (with no noise) presents also some new effects re- 
spect to the corresponding homogeneous version. In this 
section we discuss some examples related to the influence 
of initial conditions. 



A. Cultural interaction and dispersal without 
growth 

Differently from a homogeneous model, such as the AS 
model, described by ordinary differential equations, in 
the spatial model the evolution in time and space depends 
on the form of the initial population densities /i(r,io) 
(c.f. Sec. Ill A|l . While this is a standard mathematical 
property, it has a relevant meaning in terms of geographi- 
cal and historical conditions. We consider the discretized 
Eqs. (ini) on a square lattice with reflecting boundary con- 
ditions, assuming F = (no advection) and a — (no 
growth). Here and in the other simulations units arc cho- 
sen to have a constant diffusion coefficient D = 1. Such 
a model represents a simplified version of a region that 
is isolated and geographically homogeneous. We use a 
lattice of size 50 x 50 with Ax = Ay = l; the time step is 
Ai = 0.01 and the reaction constant fc = 2000. 

In order to check the consistency of the minimal spa- 
tial model, we have verified the AS model predictions by 
choosing uniform initial conditions, fi{x,y,tQ) = const. 
In this case the diffusion terms in Eqs. ([3]) or (IS|) are 
zero and the distributions remain uniform at any time i; 



this is the only case in which integrating over the space 
coordinates exactly gives back the AS model ([1]). Uni- 
form initial conditions for /i and /2 represent a historical 
moment when the populations 1 and 2 were broadly dis- 
tributed across the territory. Population 1 is observed 
to disappear whenever the initial density ratio /1//2 is 
lower than the critical value N^ /N2 given by condition 

©■ 

In the non-uniform case, we have found for various val- 
ues of parameters that a broader initial distribution rep- 
resents a disadvantage when population growth is negli- 
gible (a « 0). We illustrate this effect for two languages 
with status si =0.55 and S2 = l — si = 0.45 (si > S2) and 
equal initial total population sizes Ni{to) = N2{tQ) — l/2. 
Since q; = 0, the total number of speakers is conserved and 
for simplicity we normalize it to one, Ni(t) + N2(t) = 1. 
Thus Ni (t) represents here the fraction of speakers of the 
i-th language at time t. With such parameters and ini- 
tial conditions, language 1 would be clearly favored in 
the uniform case (or in the AS model) while language 
2 would disappear. However, this does not necessarily 
happen when space dimensions are taken into account. 

Let us investigate the situation when populations 1 and 
2 are initially distributed according to Gaussian densities. 



Mx,y,to) 



2TTaf 



(x-Xi)"^ + [y-Vif 



2oi 



(7) 



with i = 1,2, Xi — yi — 25, i.e., the average positions are 
located in the center of the simulation area (see Fig. [1] 
top). We remark that, given the symmetry of the ini- 
tial configuration, using reflecting or periodic boundary 
conditions leads to perfectly equivalent results. For pop- 
ulation 1 we assume tJi = 10 in both examples A and B, 
whereas for population 2 we assign cr2 = 1-75 in example 
A and cr2 = 3 in example B. The particular initial con- 
figurations assumed can be interpreted from a historical 
point of view as the sudden appearance of a high pop- 
ulation density of speakers 2 in the center (of mass) of 
population 1 2]|. One is then interested in predicting 
the final configuration, i.e., which language will even- 
tually prevail. Surprisingly, even when the two initial 
population sizes are equal, it is not the language with a 
higher status which necessarily survives. This situation 
is illustrated by example A in Fig.[T] In this case, the ini- 
tial population distribution, with (72 = 1-75, is sufffciently 
narrow and the associated high population density fa- 
vors the 1^2 reaction, see the function R{fi,f2) given 
by Eq. (|3|). Eventually language 1 disappears despite 
its higher status and the same initial population. Start- 
ing from a wider initial density of population 2, tT2 = 3, 
while all other parameters maintain the same values as in 
example A, the opposite final configuration is recovered, 
i.e., it is population 2 which now disappears (see example 
B in Fig. [J). 

In Fig. [2] we compare the total populations Ni{t) = 
l — N2{t) for the same examples A and B of Fig. [T] As 
one can notice, the total population TVi of example B be- 
comes smaller than N2 immediately at t > ^o = and 
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FIG. 3: Comparison of the evolution of population densities fi{x,y,t) and f2{x,y,t) (columns 1 and 2) for two languages 
with status si = l — S2 = 0.55 for different widths of the initial distribution fi{x,y,to). All other parameters are the same, see 
text. Example C: initial distribution fi{x,y,to) localized within a radius Ri =15. Example D: uniform initial distribution 
/i(a;,y,io)=const. 



remains smaller until t « 400; however, eventually pop- 
ulation 1 prevails. This is possible since language 1 has 
a higher status, si > S2- The snapshot at time i~ 100 
in Fig. [T] shows that population densities have become 
almost uniform by that time, so that one can use the AS 
model to estimate the critical ratio for the survival of 
language 1. From Eq. Q, one obtains A^j^/iVI « 0.5122 
for si — I — S2 = 0.55, corresponding to a critical fraction 
iVf « 0.3387. In Fig. [21 in the example B, the surviving 
population Ni{t) > N* at any time t , while in the ex- 
ample A Ni{t) crosses the line A^i = N^, beginning its 
irreversible decrease. 



B. Cultural interaction and dispersal with growth 

In this subsection we show that taking into account 
population growth, for sufficiently large values of a, a 
larger initial spreading can lead to the survival of the lan- 
guage, on the contrary to the situation with no growth, 



considered above. In Fig. [3] we present two illustrative 
examples, C and D, which differ from each other with 
respect to the initial spread of population 1 . 

Performing the numerical simulations, we have used a 
100x100 simulation area with Ax — Ay — 0.5 and periodic 
boundary conditions. The time step is Ai = 0.001 and the 
reaction constant k — 1000. Regarding the growth term, 
a rate a = 0.03 and a carrying capacity i^ = 0.1 have been 
used. The initial densities have been chosen of the form 



fi{x,y,h) = 



M,. 



2 Mi 



ri{x,y) < Ri, 



l+oxp{-[r.(:r,y)-fl.]V2af} ' rt{x,y)>R^, 



where ri{x,y) — ^J{x ■ 



(8) 

Y ~^ {v ~ ViY is the distance 
between position {x,y) and the average position {xi,yi) 
of the population density fi. The average positions 
have been chosen in the middle of the simulation area, 
^i = Hi = 25. The function ([5]) defines a population 
localized around {xi,yi) within a radius Ri with con- 
stant density fi = A/i, smoothly decreasing to zero for 





FIG. 4: Time evolution of the total populations Ni{t) and 
N2{t) corresponding to the examples C (left) and D (right) of 
Fig. El 



ri{x,y) > Ri on a scale <Ti. We have ensured that the 
density ([8]) is nowhere larger than the carrying capac- 
ityj fi{x,y,to) < K. In both examples C and D, pop- 
ulation 2 starts from the same initial density ([8]) with 
Af2 = 0.0163647, R2 = 4, and 0-2 = 2. Instead, for popu- 
lation 1 two different initial conditions are used. In ex- 
ample C it is given by the density ([5]), with parameters 
Afi = 0.0020048, i?i = 10, and a^ = 3; In example D, the 
parameters are M = 0.000548599, i?i = 25 and cti = 10; 
with such large values of i?i and ai , the initial density /i 
in example D is practically uniform. The parameters Mi 
and Af2 have been chosen in order to have the same total 
initial populations in both examples C and D, namely, 
Ni{to) = 1.36 and A^2(io) = 2.64. 

From Fig. [3] one observes that, differently from what 
happens in examples A and B (Fig. [T]) , increasing the ini- 
tial spread of a population density favors its survival: a 
more localized initial density /i (example C in Fig. [3]) is 
observed to disappear, while a uniform initial condition 
for /i (example D in Fig. [3]) leads to its survival. This 
can be traced back to the fact that in example D a higher 
growth of population 1 takes place in the peripheral re- 
gions, where the population density /2 is negligible. The 
time evolution of the total population Ni{t) and N2{t) 
for examples C and D are depicted in Fig. 21 



IV. INFLUENCE OF BOUNDARY 
CONDITIONS 

As shown for instance in the study of the three-state 
voter model in Ref. |ll| , boundary conditions can have a 
crucial influence on the competition process between cul- 
tural traits. Here we investigate the problem of language 
competition comparing the influence of reflecting and pe- 
riodic boundary conditions. Using reflecting boundaries 
one can model a geographical area which is isolated, i.e., 
it is not possible to enter or leave it. Periodic boundaries 
are a numerically convenient way to simulate an open re- 
gion in the middle of a much larger accessible area. We 
found that if the growth rate is negligible (a « 0) the 
vicinity of reflecting boundary conditions favors the sur- 
vival of a language. For high growth rates, the effect of 
different boundaries is less appreciable. Here we present 
the results for the case with no growth (a = 0) . In Fig. [5] 
we compare the time evolutions of population densities 
fi(x,y,t) for two languages in the presence of reflecting 
and periodic boundaries. For the language status, the 
values si = 0.55 and S2 = 0.45 have been assigned. In 
both examples the initial distributions fi{x, y, to) are as- 
sumed to have the Gaussian shape defined by Eqs. ([7]), 
with xi — 20, j/i — 30, <Ji = 3, for population 1, and 
X2 = 45, 2/2 = 5, (72 = 1, for population 2; the normaliza- 
tion constants are A^i(to) = 0.37 and N2{to) = 0.63. The 
size of the simulation area is 50 x 50 with Aa:: = Ay = 1, 
the time step is At — 0.01, and the reaction constant 
k = 1000. 

From Eq. ^ one obtains, for si = 1 — S2 = 0.55, that 
the critical fraction ensuring survival of language 1 is 
N* w 0.36 {N* = l-N*fiiOM). Thus, the initial popu- 
lation fractions A^i =0.37 and iV2 = 0.63 used here would 
give a slight advantage to population 1 in the uniform 
case. Instead, as one can see from Fig. [SJ language 1 dis- 
appears with reflecting boundary conditions (Fig.[5l left), 
while it prevails if periodic boundary conditions are used 
(Fig. m right). This effect is due to the fact that re- 
flecting boundaries bounce back a part of population 2 
located near the boundary increasing the corresponding 
density /2. With open or periodic boundaries, population 
2 would spread and its density /2 decrease. This in turn 
would lower the term representing the 2 — > 1 language 
switching rate in the function i?(/i, /2). Instead, a higher 
density /2 favors the switching of speakers 1 to language 
2. In Fig. [6] the total populations Ni{t) ^ l-iV2(i) for 
the two examples with reflecting and periodic boundary 
conditions are plotted. 



V. GEOGRAPHICAL BARRIER 

In real situations the coexistence of more than one 
langua ge i n neighboring areas for long times is often ob- 
served \V^. Some models, such as the AS model, describe 
situations in which after a relatively short time only one 
of the competing languages survives. Survival of different 
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FIG. 5: Comparison of the evolution of population densities fi{x,y,t) and f2{x,y,t) (columns 1 and 2 respectively) with 
language status si = 0.45 = l — S2 for reflecting (RBC, left part) and periodic (PBC, right part) boundary conditions, while all 
other parameters are the same. See text for details. 
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FIG. 6; Time evolution of the total population Ni{t) = 1- 
N2{t) for the examples RBC and PBC of Fig. [3 



languages in the same area was shown in Ref. [18[ to be 
possible when population growth is taken into account. 
Another mechanism was suggested in Ref. [19], based on 
the similarity of the competing languages. In the frame- 
work of the bit-string language evolution model, the in- 
fluence of a barrier on language diffusion and evolution 
was considered by Schulze and Stauffer, who showed that 
two different languages can exist on opposite sides of the 



barrier [g, Q. 

In the present paper we concentrate on the role of 
purely geographical factors. The scheme employed is dif- 
ferent from that of Ref. [S] , where the coexistence of two 
languages in neighboring regions was made possible by 
a barrier (geographical boundary or political border) af- 
fecting the form of the switching rate i?(/i, /2). Instead, 
the mechanism described in the examples presented in 
Sees. IVl and IVTl below is based on the presence of a geo- 
graphical barrier which influences solely population dis- 
persal, while the cultural interaction remains the same 
as in the AS model. 

As a first example of geographical inhomogeneity, we 
consider a barrier, representing e.g. a mountain chain, 
which divides the accessible area into two regions. For 
the sake of simplicity we model the problem in one di- 
mension. We assume that the populations of speakers of 
language 1 and 2 are initially localized on the opposite 
sides of the barrier, as depicted in Fig. [7] top. We are in- 
terested in the influence of the barrier on the asymptotic 
state. In order to answer this question, we have evolved 
the population densities fi{x,t) according to Eqs. ([3]), 
where V = d/dx; we do not take into account popula- 
tion growth (a = 0). Numerical integration is performed 
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FIG. 8: Time evolution of total populations Ni{t) and N2{t) 
in the presence of a barrier, example of Fig. [T] 



?s 





FIG. 7: Evolution of the speaker population density fi{x,t) 
(continuous line, left axis) and f2{x,t) (dashed line, left axis) 
in the presence of a barrier U{x), drawn at i = (top, gray 
area, right axis). Language 1 is favored both in status (si = 
0.6 = 1-S2) and initial population [iVi(io) =0.8 = l-iV2(io)]- 
At times i = 10 (middle) and t — 3000 (bottom) also the local 
reaction rate R{fi{x,t), f2{x,t)) is depicted (dashed-dotted 
line, right axis). In the asymptotic state (bottom) both lan- 
guages survive, being localized on the opposite sides of the 
barrier. See text for further details. 



through the Crank-Nicolson method [l5|, assuming re- 
flecting boundary conditions. A space step Ax = 0.05, 
a time step 6t = 0.001, and a reaction constant k = 200, 
are used. The barrier is modeled through the force field 
F = F{x) = —dU{x)/dx, where the potential U{x), de- 
picted in Fig. [7] top at t = to = (gray area), is 



Uix) 



Uo 



1 + exp[-(a; - Xa)/(Tu ] + exp[(a: - Xb)/cru ] 



(9) 

This function represents a step of height Uq located in 
the interval between Xa and xi, {xa < xt,) and going to 
zero outside of it on a scale length au . The parameters 



FIG. 9; Schematic map (top) and potential energy landscape 
representing the sea barrier (bottom) for the three-island con- 
figuration. Islands A, B, and C have a circular shape with the 
same radius, defined by the effective potential (|ll|l . Islands 
B and C are closer to each other than islands A and B. See 
text for further details. 



employed are: C/o = 3, Xa = 12.5, Xb — ^7.5, and cjy^O.l. 
For the initial densities fi{x,to) (Fig. [7] top, left side of 
the barrier, continuous line) and f2{x,to) (Fig. [7] top, 
right side of the barrier, dashed line) a Gaussian shape 
was assumed, 



Ji[x, to) = exp[-(a; 



2n(Ti 



x^r/2af 



(10) 



We choose ai — a2 — I and for the average coordinates 
xi=5 and X2 = 25, located symmetrically respect to the 
barrier and the reflecting boundaries at x = and a; = 30. 
The status of language 1 is si = 1 — S2 = 0.6 and the ini- 
tial population fraction Ni{to) — I — A^2(io) = 0.8. Thus, 
language 1 is highly favored, regarding both status and 
initial total population. In fact, if there is no potential 
barrier, language 2 is observed to disappear, as it is easy 
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FIG. 10: Evolution of population density fi{x,y,t) (left column) and f2{x,y,t) (right column). Notice the temporary presence 
of population 1 on (the central) island B at t~60. Despite the lower status of language 2, geographical inhomogeneities favor 
its immigration to the central island B. See text for the details. 



to guess. In the presence of the potential barrier ^ the 
evolution of the population densities is shown in Fig. [T] 
both languages survive in the asymptotic limit on the op- 
posite sides of the barrier. The time evolution of the total 
populations Ni{t) is shovifn in Fig.[8l The survival of lan- 
guage 2 on the right side of the barrier can be explained 
as follows: the potential barrier modulates diffusion to- 
ward both the directions and in particular it decreases 
the flux of population 1 from the left toward the right 
region. This in turn causes the term k sifi{x,t)f2{x,t) 
in the reaction rate ([4]) to remain very small in the right 
region, so that the local 2^1 language switching rate is 
negligible respect to that of the complementary process 
1^2. 

The survival of both languages observed in the exam- 
ple discussed above can be ascribed to the interplay be- 
tween historical conditions (initial localization of the two 
speakers communities on opposite sides of the barrier) 
and geographical constraints (presence of the barrier) . In 
Fig.[7]center and bottom also the switching rate R(fi, /2) 
is depicted with a dash-dotted line. The largest values of 
\R\ are located close to the barrier borders, where speak- 
ers coming from the other side meet the local speakers 
and switch to the local language. In Fig. [8] one can also 
notice that the asymptotic total population values are 
equal -/Vi(i— s-oo) = N2{t^oo) — 1/2, as a consequence 
of the identical dispersal properties assumed for the two 
populations and the symmetrical geometry of the system. 



VI. IMMIGRATION TO AN ISLAND 

Let us now see, what happens if we consider the spread- 
ing of two languages toward the same initially empty re- 
gion. To make an example, we study two islands A and 



C, initially colonized by populations speaking language 
1 and 2, respectively. Language 1 has a higher status 
si = 1 — S2 = 0.6. Between the two islands A and C is 
located a third island B, which is empty. For simplicity, 
we choose for the islands a circular shape with identi- 
cal radius R; the centers of the islands are located on a 
line. The situation is illustrated in Fig. [9] top. The sys- 
tem has been studied on a rectangular simulation area 
of sides Lx = 35 and Ly = 11. The population densities 
were evolved through the explicit Euler algorithm ^ on 
a 350 X 110 lattice with steps Sx = Sy = 0.1, using a time 
step 6t = 0.001 and a rate constant k — 2000. 

Speakers can freely diffuse inside the islands, but have 
to overcome a barrier in order to cross the sea and reach 
other islands. The effective potential U{x,y) modeling 
the barrier due to the sea also defines the island shapes; 
it is depicted in the lower part of Fig. [9] and is given by 

U(x v)^( J7oexp{-[rj(x,y)-i?]V2cr2,}, r^ < R, 
^ ' ^'^ \ t/o , otherwise. 

(11) 
Here j = A,B,C labels the islands and rj = rj{x,y) — 



V^ 



^j)' 



{y — yjY is the distance between the 



generic position {x,y) and the center {xj,yj) of the j- 
th island; i.e., a value rj{x,y) < R for a given j corre- 
sponds to a point {x,y) inside the j'-th island, while if 
rj{x,y) > R for all j = A, B,C the point {x,y) is on 
the sea. The location of the three islands can be rec- 
ognized in Fig. [9] bottom as the three potential wells, 
while the sea is represented by the plateau U{x,y) = Uo- 
The parameter values used are: au ~ 2, defining the 
smoothness of the potential step, Uq — 4, for the bar- 
rier height, and i? = 5, for the island radius; the coor- 
dinates of the centers of the island are xa = 6, xb = 19, 
xq = 29.25, and J/a = 2/b = yc = 6. The initial total 
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FIG. 11: Time evolution of the speaker community 1 (con- 
tinuous line) and 2 (dotted line) for the three island example. 
Top: total populations. Bottom: partial population 1 on is- 
land A and B and partial population 2 on island B and C 
(the partial populations not shown are negligible). Notice 
that both populations 1 and 2 are present on the central is- 
land B until t « 100. 



populations on islands A and C have been assigned the 
same value iVi(to) = .^2(^0) = 1/2- The population den- 
sities fi{x,y,tQ) have the same Gaussian shape (O with 
CTi = 172 = 2 and average coordinates Xi and j/i coinciding 
with the center coordinates of the respective islands. 

Language 1, with the higher status si = 0.55, would be 
favored in the absence of barriers. We would like to know 
if on island B eventually prevails language 1 or 2. It is 
clear that if the central island B is located symmetrically 
between islands A and C then it is language 1 that in 
the end will prevail (also on island A) due to its higher 
status. Population 2 may still survive on island C, where 
it was initially dominating, if the sea represents a large 
enough barrier, due to the effect described in the previous 
section which would transform islands C into a refugium. 
If not, language 1 may prevail finally also on island C. 

On the other hand, if it is easier to cross the B-C rather 
than the A-B channel, e.g. if island B is closer to island 
C than to A, then language 2 may spread first on island 
B and then maintain its superiority thanks to the barrier 
due to the presence of the A-B channel. We have studied 



the problem for various values of the distance dec be- 
tween island B and C, while keeping constant the other 
distance cZab- The example illustrated in Figs. [TUl and 
111! corresponds to a value cZab = 3 and to a much smaller 
distance dec = 0.25 between islands B and C. Figure [TOl 
shows how populations 1 and 2 disperse over the neigh- 
boring island starting from their initial locations. Due to 
the geographical asymmetry favoring dispersal of popu- 
lation 2, there is a much larger flow of population 2 from 
island B to C than population 1 from island A to B. This 
in turn leads to a rapid spreading of language 2 on island 
B. Once language 2, with a lower status, dominates on 
island B, the wide sea barrier between island A and B 
will maintain the advantage gained by language 2. Fig- 
ures [TU] and [11] also show a small presence of language 1 
on the central island B limited to a short time interval. 



VII. CONCLUSION 

Dispersal in space and time of two languages or cul- 
tural traits competing in the same geographical area has 
been studied through an extended language competition 
model based on the one proposed by Abrams and Stro- 
gatz. We have discussed how initial and boundary con- 
ditions, as well as geographical inhomogeneities, have a 
relevant (even drastic) meaning for language spreading 
and competition. 

We have observed various examples where a language, 
which in the corresponding homogeneous model would 
disappear, actually survives. Namely, in a homogeneous 
model (without space dimensions), for given values of the 
parameters, the evolution of a population is determined 
by the total initial population. In the diffusion model 
studied here, the same total population of speakers, ini- 
tially distributed in space in different ways, may evolve 
toward opposite asymptotic scenarios. Another result of 
our investigation shows that, when growth is negligible, 
a language, whose dispersal is more affected by the geo- 
graphical boundaries, is favored respect to the case when 
no limiting boundaries are present. In the presence of 
geographical inhomogeneities, modeled as potential en- 
ergy barriers, languages can survive in different regions, 
despite the possibly lower status and smaller initial pop- 
ulations. The effects discussed in the present paper are 
purely geographical, in the sense that they are related to 
the diffusion processes and to the modulation of diffusion 
due to inhomogeneities of the background. They influ- 
ence culture spreading only indirectly and are not due to 
a change in the cultural interaction law, differently from 
the model introduced in Ref. Q. 

The diffusion model and the highly idealized examples 
presented in this paper are intended as a first step toward 
a quantitative description of the space-time diffusion of 
language and cultural traits. Our study will hopefully be 
useful in solving some of the many challenging problems 
regarding language diversity [17]. At the same time, a 
more detailed understanding of the mechanisms underly- 
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ing culture transmission could be valuable concerning the 
alarming rate of disappearance of cultural and linguistic 
diversity. 
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